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The sliding columnar phase is a new liquid-crystalline phase of matter composed of two-dimensional 
smectic lattices stacked one on top of the other. This phase is characterized by strong orientational 
but weak positional correlations between lattices in neighboring layers and a vanishing shear modulus 
for sliding lattices relative to each other. A simplified elasticity theory of the phase only allows 
intralayer fiuctuations of the columns and has three important elastic constants: the compression, 
rotation, and bending moduli, B, Ky, and K. The rotationally invariant theory contains anharmonic 
terms that lead to long wavelength renormalizations of the elastic constants similar to the Grinstein- 
Pelcovits renormalization of the elastic constants in smectic liquid crystals. We calculate these 
renormalizations at the critical dimension d = 3 and find that Ky{q) ~ K^^'^{q) ~ B~^^^(q) ~ 
(ln(l/g))^'"', where g is a wavenumber. The behavior of B, Ky, and K in a. model that includes 
fluctuations perpendicular to the layers is identical to that of the simple model with rigid layers. We 
use dimensional regularization rather than a hard-cutoff renormalization scheme because ambiguities 
arise in the one-loop integrals with a finite cutoff. 
Pacs; 61.30.Cz, 64.60.Ak, 87.10.+e. 



o 



> 

00 

in 
o 

00 

a\ 



I 

O 

o 



X 



I. INTRODUCTION 

DNA, which is a semi-flexible polymer, and cationic 
lipids in solution form complexes in which the negative 
charge of the DNA is nearly compensated by the positive 
charge of the lipids. These complexes are under intensive 
study as possible nonviral carriers of DNA to cell nuclei 
for gene therapy j^. Radler, et al. have shown that under 
appropriate conditions the complexes self-assemble into 
multi-lamellar structures ||]. The lipids form stacked 
bilayer sheets with DNA molecules intercalated in the 
galleries between the bilayers as shown in Fig. ^ Each 
gallery is thick enough to accommodate only one DNA 
molecule and its hydration layer. Within each gallery, 
DNA molecules adopt a linear rather than a coiled con- 
figuration and form a regularly spaced parallel array that 
in the absence of couplings to DNA in other galleries is 
a two-dimensional smectic liquid crystal Q . The experi- 
mentally determined X-ray structure factor of these com- 
plexes is well modeled by a stack of weakly coupled 2D 
smectic lattices [||. 

Two recent theoretical papers |^,^ have pointed out 
that weakly coupled 2D smectic lattices form a new phase 
of matter, the sliding columnar phase. This phase is char- 
acterized by strong orientational correlations but weak 
positional correlations between smectic lattices. All lat- 
tices are aligned on average along a common direction 
(the a;-axis in Fig. |^), but their relative positions decorre- 
late exponentially with distance between smectic lattices. 
With sufficiently strong coupling between galleries, long- 
range positional correlations between smectic layers de- 
velop, and the system becomes an anisotropic columnar 
phase with a two-dimensional DNA lattice in the plane 
perpendicular to the direction of DNA alignment. The 



sliding columnar phase, on the other hand, is what the 
columnar phase becomes when coupling between galleries 
becomes so weak that DNA lattices can slide freely across 
each other. It has no shear modulus resisting relative dis- 
placements of DNA lattices, but it does have a rotation 
modulus resisting their relative rotation. Dislocations 
may melt the sliding columnar phase to an anisotropic 
nematic lamellar phase at length scales longer than an 
in-plane dislocation unbinding length ||^. It is possible, 
however, to choose interlayer interactions so that the slid- 
ing columnar phase is the stable equilibrium phase at all 
length scales B. 




FIG. I. Picture of the idealized sliding columnar phase. 
The DNA columns are sandwiched between planar lipid bi- 
layer sheets. The bilayer planes are stacked in the j/-direction 
with spacing a. The DNA columns are oriented in the x di- 
rection, and, within each layer, the columns are separated by 
d. The positions of columns in neighboring layers are uncor- 
related. 
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This paper will investigate the nonlinear elasticity of 
the sliding columnar (SC) phase. Its principal purpose 
is to show that the nonlinear strains lead to a Grinstein- 
Pclcovits renormalization of the elastic constants ^ and 
not, as one could imagine, to the destruction of the slid- 
ing columnar phase itself. The lipid bilayers, which we 
take to be aligned on average parallel to the xz plane 
as shown in Fig. |l], fluctuate like bilayers in any lamellar 
phase. To understand correlations and fluctuations of the 
DNA smectic lattices, it is convenient to consider first a 
model in which the lipid bilayers are rigid planes with 
no fluctuations in the y-direction. In this case, displace- 
ments of the DNA lattices, which are aligned on average 
along the x-direction, are restricted to the z-direction. 

The rotationally invariant Landau-Ginzburg- Wilson 
Hamiltonian in units of fc^T for this system is 

n=^ld'x [Bui + Ky{d^dyu.f + K{dlu.f] , (1) 

where Ky, and K are the compression, rotation, and 
bending moduli divided by fc^T and 



(2) 



is the nonlinear Eulerian strain appropriate for the sliding 
columnar phase. Note that Ti. is invariant under 



w'^(x) ^ M^(x) + f{y). 



(3) 



It is this fact that ensures that nonlinearities do not de- 
stroy the sliding columnar phase. 

The rotationally invariant strain Uzz introduces an- 
harmonic terms into the Hamiltonian that lead to a 
Grinstein-Pelcovits renormalization of B, Ky, and K. 
The renormalized moduli scale logarithmically with q at 
long wavelengths: 



Kyiq) ~ K'/\q) ~ B-'/^q) 



ln(^ 



1/4 



(4) 



where /i is a large momentum cutoff. A complete model 
for the sliding columnar phase allows both lipid bilayers 
and smectic lattices to fluctuate. This model also ex- 
hibits Grinstein-Pelcovits renormalization of the elastic 
constants. Table | lists the exponents describing singu- 
larities in the elastic constants for both the 3D smectic 
and sliding columnar phases. 

TABLE I. Comparison of the logarithmic scaling expo- 
nents for the elastic moduli of the 3D smectic and sliding 
columnar phases. At long wavelengths the elastic moduli for 
both phases scale as ln"[l/g] with a given below. 



Phase 


B 


K 


Ky 


3D smectic 


-4/5 


2/5 




sliding columnar 


-3/4 


1/2 


1/4 



The evaluation of the above renormalization presented 
some unexpected difficulties. The continuum Hamilto- 
nian in is formally invariant under arbitrary global 
rotations. However, the introduction of a hard cutoff 
breaks this rotational invariance just as the introduc- 
tion of a similar cutoff breaks gauge invariance in gauge 
Hamiltonians [Q. Nevertheless, hard cutoff RG proce- 
dures can with care be applied successfully to Hamilto- 
nians with gauge or rotation symmetries In- 
deed, the original Grinstein-Pelcovits calculation of the 
logarithmic renormalization of the smectic-A elastic con- 
stants used a hard-cutoff |^ . When we applied the popu- 
lar momentum-shell hard cutoff RG procedure to the 
nonlinearities in the sliding columnar phase, we encoun- 
tered ambiguities that we were unable to resolve. We 
found that the values of the one-loop diagrams depended 
on whether the external momentum was added to the top 
or the bottom part of the internal loop. Similar difficul- 
ties are not encountered in the Grinstein-Pelcovits cal- 
culation. To eliminate these ambiguities, we switched to 
the dimensional regularization procedure which explic- 
itly preserves rotational invariance because the cutoffs 
are infinite |lj] . 

The remainder of the paper will be organized as fol- 
lows: we first rederive the results of Grinstein and Pel- 
covits in Sec. O using dimensional regularization. Then 



in Sec. Ill, we calculate the renormalization of the sliding 
columnar elastic constants of the simplified theory using 
the same scheme. In Sec. [V we relax the constraint of 



rigid membranes and show that the membrane fluctua- 
tions do not modify the scaling behavior of the elastic 
moduli of the rigid theory. In Appendices ^ and we 
evaluate the one-loop diagrams for the 3D smectic and 
simplified sliding columnar theories. In Appendix ^ we 
show that ambiguities arise when a finite cutoff is im- 
plemented to calculate the loop diagrams of the sliding 
columnar theory. Finally, in Appendix ^ we derive the 
nonlinear strains required for the rotationally invariant 
theory of the sliding columnar phase in the presence of 
fluctuating membranes. 



II. RG ANALYSIS OF THE 3D SMECTIC 

The rotationally invariant elasticity theory for a smec- 
tic liquid crystal contains nonlinear terms that renor- 
malize the elastic constants of the harmonic theory for 
all dimensions below three. Grinstein and Pelcovits cal- 
culated the corrections to the elastic constants of a 3D 
smectic using an RG analysis with a finite wavenumber 
cutoff 1^. They found that the corrections to both the 
compression and bending moduli are logarithmic in the 
wavenumber q with the former scaling to zero and the lat- 
ter scaling to inflnity at long wavelengths. Application 
of a hard cutoff RG procedure to the sliding columnar 
phase leads to ambiguities with no obvious resolution. 
(See Appendix O.) We, therefore, employ a dimensional 
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regularization procedure that sends the cutoff to infin- 
ity and thereby preserves rotational invariance. In this 
section we rederive the Grinstein-Pelcovits resuhs for a 
3D smectic using dimensional regularization. This estab- 
lishes the language needed to calculate the renormaliza- 
tion in the sliding columnar phase. 



A. Rotationally Invariant Theory 

A smectic in d dimensions is characterized by a mass- 
density wave with period P — 2ti /q^ along one dimension 
and by fluid- like order in the other d—1 dimensions. The 
phase of the mass density wave at point x = (xj^ , z) is 
qo'"(x). The elastic Landau-Ginzburg- Wilson Hamilto- 
nian for a smectic in units of fc^T is 



1 



(5) 



where V is the gradient operator in the d—1 subspace 
spanned by xx and Bsm and ifsm are, respectively, the 
compression and bending moduli divided by fc^T. The 
nonlinear Eulerian strain Uzz = dzU — (l/2)(Vu)'^ is in- 
variant with respect to uniform, rigid rotations of the 
smectic layers. Below we will drop the (dzu)'^ term in Uzz 
since its inclusion leads to nonlinear terms that are irrel- 
evant in the RG sense with respect to the two quadratic 
terms in (|). Therefore, we will take 



(6) 



B. Engineering Dimensions 

To implement our RG procedure it is convenient to 
rescale parameters so that Bgm is replaced by unity and 
the nonlinear form of Uzz is preserved. To this end, we 
scale u and x as follows: 



u — LuU, z — LzZ, and xx = xj 



(7) 



Note that x_l does not rescale. Under these rescalings we 
obtain 



1 



(8) 



We require Uzz = Auzz where Uzz = d^u — (l/2)(V~u)^ 
is the rescaled nonlinear strain with the same form as 
(|). This yields ^ L-^ and A = Ll. The coefficient 
of in the rescaled Hamiltonian is set to one with the 
choice 



T — R"l/3 

The rescaled theory then becomes 



(9) 



n 



d'^'a 



with 



B, 



1/3 



(10) 



(11) 



For the remainder of Sec. || we will use the free energy 
in but drop the tilde on the scaled variables. 

We determine the dimensions of the scaled variables 
using the engineering dimensions of Bsm and Ksm- The 
dimension dA determines how A scales with length L: 
[A\ = L'^^ . From the respective dimensions ds^^ = ~d 
and — 2 — d of Bsm and Ksm, we obtain [Lu] — 

[L-i] = L''/^. Using these we find the following for the 
dimensions of the scaled variables and the parameter w: 



Lu 
L 



- L, and [w] 



L 

Tz 



= 1^+"/^ (12) 



where e — i — d. Using these definitions one can eas- 
ily verify that both terms in din) are dimensionless. \w\ 
scales as ^'^'^/^ where [^] = i , and it is, therefore, a 
relevant variable below d — i. We introduce a dimen- 
sionless coupling constant go via 



w 



(foA*) 



2e/3 



(13) 



to display explicitly the length dependence of w. The 
dimensions of the coefficients of the {dzu)^ and [dzuY 
terms we omitted are, respectively, 2c?/3 and Ad/i. These 
terms are irrelevant and will be ignored in what follows. 

The engineering dimensions in ( p^ ) imply that there is 
an invariance of TL under the transformation /i — > /i6 and 



u(x_L; z) — 6'*"u'(x^, z'), 
where x'^ = 6^^xj^ and z' = b^^^^*^/^^ z, i.t 



(14) 



(15) 



This in turn implies a scaling form for the position cor- 
relation function G(xj^,z) = (u(x^, z)u(0, 0)) and its 
Fourier transform G(q). We find 



G(x^, z, w) = 52(i-'^/3)G(x'^, z', wh"'/^), 



(16) 



and from this we obtain the vertex hmction r(q) = 

r(q^,g„H =?>-'^'+'/'^r(6q^,6Wg.,z«62'/^). (17) 
When = 3 this reduces to the scaling form 



r(q_L,gz,w) = q\T [l, 

which the harmonic vertex function F = 
isfies. 



(18) 



w sat- 
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C. RG Procedure 

To calculate renormalized quantities, we seek a mul- 
tiplicative procedure that yields a renormalized Hamil- 
tonian with the same form as the original Hamiltonian, 
i.e., a Hamiltonian that is a function of a renormalized 
nonlinear strain with the same form as @. To preserve 
the form of the strains, it is necessary to rescale fields 
and lengths simultaneously. The rescaling that produced 
( p^ ) shows that the form of Uzz is preserved if the rescal- 
ing coefficients of u and z are inverses of each other. We, 
therefore, introduce a renormalization constant Z and a 
renormalized displacement u' such that 



(19) 



This imphcs that u^zix) = Z^^^u'^zi^')- We also in- 
troduce a unitless renormalized coupling constant g and 



renormalization constant Zg via 



(20) 



where /i is an arbitrary wavenumber scale. The renor- 
malized Hamiltonian then becomes 



z{u'zzY + {9^i'Zg 



-2/3 



(21) 



We now follow standard procedures to evaluate Z{g) and 
Zg{g) [Q. The renormalized Hamiltonian in ( pl| ) deter- 
mines the vertex function 



r(q) = ql + (gff 
+ (.9M' 



-2/3 4 



-2/3 



^ -2/3 



{Z - l)ql (22) 



to one-loop order, where S(q) is the one- loop diagram- 
matic contribution to r(q). We next impose the following 
conditions on the vertex function to enforce the correct 
scaling behavior: 



dV 

dql 
dT 



1 



(jj=/j2,(jj_=0 



N-2/3 



(23a) 
(23b) 



In Appendix ^ we show that the diagrammatic contri- 
butions are the following: 



dS(q) 



dQz 
dE(q) 



dqi 



167re 
(5M' 



2/3 g 

327re' 



(24a) 
(24b) 



Using the conditions on the vertex function we obtain 
the relations for the renormalization constants in terms 



of the one-loop diagrammatic corrections. The following 
relations are correct to lowest order in e: 



Z = 1 



9 



3g 

647re' 



(25a) 
(25b) 



1. Callan-Symanzik Equation 

The renormalized vertex function rr(q) satisfies a 
Callan-Symanzik (CS) equation under a change of length 
scale fj,. We obtain the renormalized elastic moduli from 
the solution to this equation. The original theory in ( [l0| ) 
did not depend on the length scale /i. We can therefore 
write the bare vertex function T in terms of the renor- 
malized vertex function and find the differential equa- 
tion obeyed by F^. Since the variables u and z scale as 
u'(x) — Z^^'^u{x') and z' — Z^l'^z, the vertex function 
must scale as 

F(qi,9.,w) = Z-i/3r^(qi,Z-i/3g„g,;,). (26) 

The CS equation is determined by the condition 
fidT/dii ~ 0. Since the renormalized vertex function can 
have explicit as well as implicit /i dependence through 
the functions Z and g, the CS equation for F^ has three 
terms: 



_d_ 

dqz 



d_ 

da. 



where 



dg_ 

dfi ' 



via) = Pig) 



dilnZ) 
dg 



Tr = 0, (27) 

(28a) 
(28b) 



and qzd/dqz = q'zd/dq'^ with q'^ = Z~^^^qz. This equa- 
tion can be integrated to yield an equation for F^ as a 
function of the length scale /i. 



rr(q_L,g2,ff,M) = exp 



r]dl' 



xTr q_L,exp 



3 7o 



rjdl' 



where /^//io 



fid/dfi 



d/dl, and 

dgil) 
dl ' 



(29) 

qz,gil),fio^ , 
(30) 



At / = we have set Tril 0) = Fr(q_L, (/z, ffo, Mo)- Now 
we must solve for P and rj in terms of g in order to obtain 
the renormalized vertex function. To find /3(g), we note 
that 



4 



dl 



dl ( 



0. 



(31) 



From this relation we find /^(g) = —e/{d{\nQ)/dg) where 
Q = gZgZ^I"^ . We can then plug in the relations for Z 
and -Eg, and we determine /3 and 77 to be the following: 



9-^9 



Vi9) = -T^f- 
Idtt 



(32a) 
(32b) 



In three dimensions e = 0. In this case, integration of 
dg/dl yields 



5(0 = 



90 



l + 5.goV(647r)' 



(33) 



where go = g{0) — w^^^. The remaining task is simple; 
we must evaluate the arguments of the exponentials in 
( p9| ) to obtain the I dependence of F^. Since g ^ l/l, the 
integral of i] will scale as In I and the exponentials of the 
integral of rj will give power-law dependence on I. We 
find that 



exp 



7]{l')dl' 



1 



— / 

647r 

- [9/9of'' 



-4/15 



(34) 



2. Renormalized Elastic Constants 

The scaling relations in ( p^ and ( p^ ) imply that F^ 
satisfies 

Triq±,q.,g,fi)^b-^[g/gof (35) 
xTr{bq±,b^ [g/gof^^^q,.,g,fiob). 

We now choose the reference length scale b — fj,Q 

2 , .4 \-l /4 _ r7_/-.M-l 



h{q)] . This implies that 



I = In 



hiq) 



(36) 



since /i//io — e'. We find the scaling form of the renor- 
malized vertex function. 



r"'[9l9.f'\i + [9/9of'ql 



(37) 



by squaring the term in the second slot of the renormal- 
ized vertex function and adding it to g~^l^ times the 
fourth power of the term in the first slot. We then plug 
in (^ ) for g and transform back to variables with dimen- 
sion to find the following expression for the renormalized 
vertex function: 



r.(q) = Bs: 



1 + ^ln 
647r 



5go 
647r 



In 



/i(q) 



-4/5 



(38) 



/j(q) 



2/5 



„l/2 -3/2 



L-sm/^sm- is a wavenum- 



and /i(q) 



where go 

(q2 + A2g4 )l/4 with A2 = 

ber A ~ 1/a associated with the short distance scale a. 
We identify the renormalized compression and bending 
moduli -Bsin(q) and -ft'sm(q) as the coefficients of the 
and terms respectively. The renormalized elastic con- 
stants scale as powers of logarithms at long wavelengths: 



i^.„,(q) ^ S-„y2(q) 



In = 



h{q) 



2/5 



(39) 



where the long wavelength regime is defined by wavenum- 
bers q that satisfy h{q) <^ A^/^ exp [— 647r/(5go)]- We see 
that Ksmiq) scales to infinity and i?sm(q) scales to zero 
as q 0. 



III. THE SLIDING COLUMNAR PHASE WITH 
RIGID LAYERS 

In this section we calculate the logarithmic corrections 
to the elastic constants for the sliding columnar phase 
using the dimensional regularization scheme employed in 
the previous section. The steps we follow for the di- 
mensional regularization of the SC phase closely resem- 
ble those followed for the dimensional regularization of 
the 3D smectic phase since the two Hamiltonians have 
similar forms. In this section we assume that each 2D 
lattice of columns is flat and only allowed to fluctuate in 
the z-direction. We relax this assumption in Sec. IV and 
find that the renormalized elastic constants are identical 
to those of the flat theory to lowest order in the coupling 
between strains in the y- and z-directions. 



A. Rotationally Invariant Theory 

The rotationally invariant elasticity theory describing 
the sliding columnar phase was derived previously in 
[^,|| . We found that a phase with weak positional correla- 
tions but strong orientational correlations between neigh- 
boring 2D smectic lattices was possible for sufficiently low 
temperatures. The strong orientational correlations re- 
quire a rotation modulus in the Landau-Ginzburg- Wilson 
Hamiltonian that assesses an energy cost for relative ro- 
tations of the lattices in addition to the compression and 
bending energy costs for a single lattice of columns. The 
Hamiltonian for the idealized sliding columnar phase in 
three dimensions and in units of ksT is 



n 



I J d'x [Bi 



+ K{dlu,f + Ky{dyd.,u,f] , (40) 
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where B, Ky, and K are the compression, rotation, 
and bending moduh divided by ksT. Symmetry per- 
mits additional terms in the Hamiltonian proportional 
to Kzy{dzdyUz)^ and KzxidzdxUzY ■ The Kzy term mea- 
sures the energy cost associated with variation in the 
DNA lattice spacing from layer to layer, and the Kzx 
term measures the energy cost associated with the vari- 
ation in the orientation with strand number of DNA 
strands within a layer. These terms are, however, sub- 
dominant to those kept in (^), and the couplings Kzy 
and Kzx are irrelevant. We will ignore them in what 
follows. The nonlinear strain Uzz is identical to the 
nonlinear strain for one layer of columns Uzz — dzUz — 
{l/2)[{dxUz)'^ + [dzUz f]. Below we will drop the {dzUz f 
term from the nonlinear strain since it leads to terms in 
the nonlinear theory that are also irrelevant with respect 
to the three harmonic terms in (^o|). Therefore, we use 
the approximate expression, 



The rescaled Hamiltonian becomes 



dzUz 



\ {dxUzf . 



(41) 



We note that Uzz and H do not possess a shear strain 
term (dyUz)^ because neighboring layers of columns can 
slide relative to one another without energy cost. The ab- 
sence of the shear energy cost is a unique feature of the 
sliding columnar elasticity theory. Because the Hamilto- 
nian lacks terms with y-derivatives alone, it is invariant 
with respect to shifts in Uz that are only a function of y. 
Hence, HlUz] — Hluz] with 



f{y)- 



(42) 



This invariance restates that there is no energy cost for 
sliding neighboring layers of columns relative to one an- 
other by an arbitrary amount. 



B. Engineering Dimensions 

We simplify the sliding columnar theory in (^0|) by 
rescaling the lengths so that B and Ky are replaced by 
unity and the nonlinear form of Uzz is preserved. We 
accomplish this by scaling Uz, y, and z but not x. To 
implement a dimensional regularization scheme it is nec- 
essary to let X become a d — 2 dimensional displacement 
in the space perpendicular to y and z. Rescaled variables 
are defined via 



Uz — Ly^Uz^ X — X, 

y — Lyij, and z — LzZ. 



(43) 



We first set L„ — to preserve the form of Uzz under 
(^). We then set the coefficients of and {dydxUz)'^ 
to unity by choosing 




and Lz = {KyBY'"^. 



(44) 



^-2 



ulz + {didyUzf +W-\dluzf 



with 



£1/2 



W = 



KK 



1/2 



(45) 

(46) 
I and 



and d = 3 — e. In the rest of this section we use 
drop the tildes. 

We determine the dimension of the scaled variables 
from the dimensions of the elastic constants in (^). The 
dimensions [B] = L''^ and [Ky] = [K] = L^-"* dictate 



[uz] = L'^^-'^)/'^^ [x] = L. 

[yf=L('*-i)/2,' [z] = L('^+i)/2, and [w] ^ L'^-^ 



(47) 



Note that [w] scales as /z^ with [/i] = L^^ and is relevant 
below d — 3. The length dependence of w is extracted 
by introducing a dimensionless coupling constant go via 
w = gofi". 

The engineering dimensions in ( p7| ) imply that the 
Hamiltonian is invariant under the transformations — > 
fib and 



(48) 



with x' = b ^x, y 



^-(d-i)/2y^ and z' = 6-(^+i)/2z, 
i.e. the Hamiltonian obeys 



(49) 



This implies that there is a scaling form for the position 
correlation function G'(x) — {uz{x.)uz{0)) and the vertex 
function F = G^^. We find that r(q) obeys the following 
scaling relation: 

r(q,H- (50) 

When d = 3 this reduces to 

riq,w)=qtT{l,qy/qx,qz/ql), (51) 
which is satisfied by the SC harmonic vertex function 



-1 4 

■w 'ql. 



C. RG Procedure 



We now follow closely the RG procedure in Sec. lIC . 
We rescale the lengths and fields, ensure that the SC 
Hamiltonian has the same form as the unsealed SC 
Hamiltonian, impose boundary conditions on the vertex 
function, and determine the renormalization constants 
in terms of the one-loop diagrammatic corrections. The 
first step in the process is to rescale lengths such that 
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the renormalized SC Hamiltonian has the same form as 
(^). To preserve the form of the nonhnear strain, the 
z and u rescahngs must be inverses of one another and 
the y rescahng is arbitrary. We, therefore, introduce two 
renormahzation constants Z and Zy such that 

u,(x) = Zi/3<(x') = Zyy, Z^l^z). (52) 

This imphes that ^^^(x) = Z^/'^u^^(x') and dxdyUz{x) = 
Z^/^Zydx'dyiu'^{x'). We also define a unitless couphng 
constant g and renormahzation constant Zg by setting 



w^gfi'Z'/'ZyZ-\ 

The renormahzed Hamiltonian then becomes 
1 



(53) 



n' 



ZZyHu'z.f+Z^/^Zy{d.x,dy,u'zf 



+ {gt^^Zy)-\dl,u',f 



(54) 



We again employ standard RG procedures to calcu- 
late Z^ Zy^ and Zg. The renormahzation constants are 
fixed once we impose the following three conditions on 
the vertex function: 

dV 



dV 



= 1 



(55) 



dV 



(Note that we have dropped the primes on the rescaled 
Hamiltonian.) The vertex function to one-loop order. 



2 2 

qxqy 



rl/3 



Zy-l] qlql 



-\t + {ZZ-' - 1) ql (56) 
(5M^)-i(Z;i-l)g4 + s(q), 



is obtained from (54) by adding and subtracting + 
q^qy + {gfj,'^)~^qx and including the one- loop diagram- 
matic contributions to the vertex function, S(q). In Ap- 
pendix H we calculate the diagrammatic contributions. 



dql 
dT, 



9 



d{qlql) 
dH 
dqt 



9 



247r2e 



9 



127r2e' 



(57a) 
(57b) 
(57c) 



to lowest order in e. From these we determine the renor- 
mahzation constants to be 

9 



Z -l-^— 
_i , 9 



(58a) 
(58b) 
(58c) 



1. Callan-Symanzik Equation 

The Callan-Symanzik equation is obtained by requir- 
ing that the original theory in ( ^ ) be independent of 
the length scale fj,. To ensure this, we set ^idT/dfi = 0. 
This can be converted into a differential equation in the 
renormalized vertex function using the following scal- 
ing relation: 



r{q,w)=Z-^/''ZyTr. (qx,Zy\y,Z-^/\, 



(59) 



From the scaling relation we determine that the CS equa- 
tion has the following four terms: 



d_ 

d 



d_ 

dg 



F. = 0, 



(60) 



where ry(g) and I3{g) were defined previously in Sec. H C 1 
and riy {g) — l3{g)d{h\Zy)/dg. The solution to ( |60| ) is 



Tr (q, .g, m) = exp 



(61) 



fe,exp 



rjydV 



qy,exp 



Tjdl' 



with Tril = 0) = Fr(q,5o,Aio) and fi/fj.(, = e'. 

The coupling constant w must be independent of the 
length scale I. This condition yields a differential equa- 
tion for the dimensionless constant g whose solution is 



5(0 



90 



l + gol/iGiT^) 



(62) 



This equation in turn determines the I dependence of rj 
and r]y since they are both proportional to g. We find 
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16^2' 



(63) 



and thus these scale as l/l at long wavelengths. 



2. Renormalized Elastic Constants 



Using ( |62| ) for g{l) and the relations for rj^g) and riy{g) 
in (|6^), we obtain the scaling form of the renormalized 
vertex function: 

Triq,g,^i)^b-*[g/go]'^' X (64) 
Tr (bqx, bqy [g/goV^^^ , b^qz [9/90]^^^ , 9, Mo^) ■ 

To set the length scale, we choose 

b = ^^o' = {ql + qlql + ^-'q^T'" ^ IMq)]"' • (65) 
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It follows that 



I = In 



(66) 



since /i and I are related via /i//io = e'. We then sub- 
stitute ( |6^ ) for g and transform back to variables with 
dimension to obtain the following expression for rr(q): 



r,(q) =B[l 



+ Ky\l 



K 1 



.90 
67r2 

go 

67r2 



In 
In 
In 



/i(q) 



-3/4 



(67) 



/i(q) 



1/4 



2 2 



1/2 



4 



where 50 = B^/^ / {KKI'^), Jl = fj./{KyBy/^, and /i(q) = 
(q2 + X^qlq^ + X'^q^^/^ with = Ky/B and = K/B. 
JI^ is an upper momentum cutoff A ^ 1/a associated 
with the short distance scale a. We can now identify 
the q dependent elastic constants and determine their 
scaling as q tends to zero. At long wavelengths such that 
/i(q) <ti A^/^exp [— 67r2/go] the In term dominates, and 
we find 



Ky{q)^K'/'{q)^B-'r^{q) 



1/3, 



In U 



h{q) 



1/4 



(68) 



We see that B{q) scales to zero and K{q) and Ky{q) 
scale to infinity as g 0. Also note in Tabic | that 
the exponents of the logarithmic power-laws of B{q) and 
K{q) are different from those of i3sm(q) and Ksm{l), but 
the signs of the respective exponents are the same. 



IV. SLIDING COLUMNAR PHASE WITH 
FLUCTUATING LIPID BILAYERS 

In the preceeding section, we considered a model for 
lamellar DNA-lipid complexes in which lipid bilayers 
were treated as rigid planes and no displacements of DNA 
lattices in the y-direction were allowed. In physically re- 
alized complexes, lipid bilayers can undergo shape fluc- 
tuations and DNA lattices can undergo j/-displacements. 
We can parameterize the shape of the nth bilayer by a 
height function hn{x,z), which in the continuum limit 
becomes /i(x) = hy^g^{x,z). The y-displacement of the 
DNA lattices in the continuum limit is Uy{x.). At long 
wavelengths the displacements h{x) and Uy(x) are locked 
together. The lock-in occurs because there is an energy 
cost for translating each lattice of columns and the lipid 
bilayers by different constant amounts in the y-direction. 
(See Fig.y.) We can, therefore, describe long wavelength 
elastic distortions and fluctuations of the sliding colum- 
nar phase in terms of a Landau-Ginzburg- Wilson elastic 
Hamiltonian expressed in terms of displacements Uz and 



nb [uy, - 2 / + KM^z? (69) 

2 



+ KyidxdyU,) 



+ KlAdAuyf + KUdluyf + 2By^UyyUz 



where Uyy and Uzz are nonlinear strains. We define Tib to 
have units of fc^T, and therefore the constants appear- 
ing in this equation are the compression and bending 
moduli divided by ksT. The first three terms in ( |69| ) 
were discussed previously in Sec. [II as the Uz theory 



for the sliding columnar phase without fluctuations of 
the lipid bilayers. The next four terms are the compres- 
sion and bending energies for an anisotropic 3D smectic 
with layers parallel to the xz plane. The bending energy 
is anisotropic due to the presence of the DNA columns. 
The final term is a coupling of the nonlinear strains Uyy 
and Uzz- 

The form of the nonlinear strains depends on whether 
Eulerian or Lagrangian coordinates are used |]l5| . We find 
it convenient to use a mixed parameterization in which 
X and z are Eulerian coordinates specifying a position in 
space and y = na is a Lagrangian coordinate specifying 
the layer number. In Appendix we derive the nonlin- 
ear strains Uzz and Uyy for this mixed parameterization. 
To quadratic order in gradients of Uy and Uz, we find 



= dzUz - 



[{d.,Uyf + [dzUyf - {dyUyf 

{d^uzf + {dzUzf - [dzUyf 



(70a) 
(70b) 



Note that the nonlinear strain 



does not contain 



the shear strain term proportional to (dyUz)'^. Thus, 
layer fluctuations do not modify the essential invariance 
u'z Uz + fiy) of the sliding columnar phase to the order 
considered here . In what follows, we will truncate the 
nonlinear strains to 



lyy Ri dyUy 



Uzz ~ dzUz - -{dxUz) 



(71a) 
(71b) 



since the other nonlinear terms are irrelevant with respect 
to the sliding columnar harmonic terms in (|69|). 

The goal of this section is to calculate the Grinstein- 
Pelcovits renormalization of the eight elastic constants 
found in the theory of the sliding columnar phase with 
lipid bilayer fluctuations. Since the nonlinear strains do 
not introduce a {dyUzY term, we do not expect the bi- 
layer fluctuations to alter the renormalization of the SG 
elastic constants in the simplified theory of the previ- 
ous section to lowest order in B^^- . We will again use 
dimensional regularization to calculate the renormaliza- 
tion. The format will closely parallel the previous SG 
calculation. We first determine which of the harmonic 
terms in ( |69[ ) are relevant and drop irrelevant terms. We 
then rescale lengths and fields, ensure that the Hamilto- 
nian retains its unsealed form, impose boundary bound- 
ary conditions on the vertex function, and calculate the 
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renormalization constants. The renormalization con- 
stants then determine the scahng form of the vertex func- 
tion. 



A. Engineering Dimensions 

We begin by rescahng the lengths and t he fields in Tii,. 
In addition to the rescalings in Sec. [II B , we also rescale 
Uy according to 



(72) 



We first impose the conditions of the previous section, 
i.e. we set the coefficients of and {dxdyUzY to unity 
and ensure that both terms in the nonlinear strain u^z 
scale the same way. As an added constraint, we set the 
coefficient of u^y to unity. These conditions fix 

, — , 1 //I _^ 




(73) 



Once we plug in these scaling lengths, the rescaled 
Hamiltonian becomes 



(dyUyY + 2v{dyUy)UzZ + VlidjUyY 

V2{didiUyY + vsidjuyY^ 



(74) 



with 



w 



Vl 



V3 



(B^)i/2 

XX \ xy ) 

Ry (K"- 

XX \ xy/ 
5^(5^)1/2 ' 



V2 



xz xy 

~ ByB^ ' 



(75) 



and 



\l/2 



{B^f/^By 



(It is again necessary to let x represent a d — 2 displace- 
ment with d = 3 — e.) The dimensions of the scaled vari- 
ables and the w and v coefficients are determined using 
( [73] ) and the dimensions of the compression and bend- 
ing moduli, [B] = L^'^ and [K] — L^^"^. (Note we have 
dropped the tildes on the scaled variables in the following 
discussion.) We find 



«2 



i^, and [vs] 



«1 



L 



5-d 



(76) 



while the dimensions of Uz, y, z, and w were given pre- 
viously in (E^). Note that v does not scale with length. 



Also note that the coefficients vi, V2, and are irrelevant 
when d = 3. We drop the irrelevant terms and arrive at 
the following simplified Hamiltonian: 



1 



Hfa = - / d^x 



+ {dxdyUzf +W-\dlUzf 

+ {dyUyY + 2v{dyUy)uzz 



(77) 



B. RG Procedure 

The present R G pro ced ure wi ll be similar to those em- 
ployed in sections II C and III C , except we now have two 



coupling constants, w and v, instead of one. We will show 
that the inclusion of v does not alter the renormalization 
of the sliding columnar elastic constants to lowest order 
in V. As before, we rescale the fields and lengths and seek 
a renormalized Hamiltonian with the same form as (|7^) . 
We scale y, z, and Uz as we did previously in (|5^) and 
Uy by Z^/^ as follows: 



Uy{-K) 



- ^1/2,/ 



(x')=Zi/2y;(x,Z,2/,Zi/3z). (78) 



The rescaled Hamiltonian looks similar to (£^) with 
two additional terms due to fluctuations of the bilayers. 
We drop the primes on the variables and find 



Z^^^'ZyidxdyUzY 



(79) 



where 



{gfi^Zy)-\d^,UzY + Z-^/'ZyZ{dyUyY 
2vZy{dyUy)Uz 

vZ„ - i;Zi/2zi/3 



(80) 



and Zg was defined previously. 

Boundary conditions imposed on the vertex functions 
Tij{q) with i,j = y,z ensure that the Hamiltonian re- 
tains its original form in (^^ after rescaling. The vertex 
function is defined by Tij{q) = G^-^(q) with Gij{x) = 
{ui{x)uj{0)). The conditions imposed on Tzz are identi- 
cal to those given in (|55|); these are augmented by two 



conditions on Tyz and Tyy. 



d{qyqz) 
dV 



2v 



(81) 



yy 



dql 



= 1. 



9 




u 



FIG. 2. Schematic diagram of the additional relevant non- 
linear term dyUyidxUz)'^ generated by the sliding columnar 
theory with lipid bilayer fluctuations. The symbols x and 
y written adjacent to the dividing lines represent x and y 
derivatives of the respective fields. The Uy field is denoted by 
a dashed line while Uz is denoted by an unbroken line. 

Once we impose these conditions on the vertex func- 
tions, we solve for the Z's in terms of the one-loop dia- 
grammatic contributions Sy , where, for instance, Y,zz is 
the one-loop correction to the vertex function T^z- The 
diagrammatic corrections arise from the quadratic term 
in Uzz- generates dzUz{dxUzY , which was already 
present in the theory with Uy — Q. The coupling oiuyy to 
Uzz generates a new nonlinear term, cLuy{dxUz)^ ■ This 



term is shown schematically in Fig. g. There are six 
new one-loop diagrams in addition to the three diagrams 
of the rigid sliding columnar theory; these are shown in 
Figs. ^ and ^. The diagrams in Fig. ^ arise from contrac- 
tions of dyUy{dxUzY with itself and the diagrams in Fig. 
arise from contractions of dyUy{dxUz)'^ with dzUz{dxUz) ■ 



(a) 






(c) 



FIG. 4. The three diagrams that can be formed by con- 
tracting dzUz{dxUzY with dyUy{dxUz)^ . The only diagram 
that contributes to the renormalization of B^^ is pictured in 
(a). The diagrams pictured in (b) and (c) contribute to the 
renormalization of both K^^ and K^y. 

The one-loop diagrammatic corrections "Ezz are easy 
to calculate since the form of the propagator Gzz is un- 
changed from its form in the rigid sUding columnar the- 
ory. The form is not changed, but the compression mod- 
ulus B is renormalized by a factor ofl — v^. The one- loop 



diagrammatic corrections to T^ 
est order in e: 



are shown below to low- 



(a) 



(b) 



(c) 




y 



FIG. 3. The three diagrams that can be formed by con- 
tracting dyUy{dxUz)^ with itself. The only diagram that con- 
tributes to the renormalization of is pictured in (a). The 
diagrams in (b) and (c) contribute to the renormalization of 



both Kxx and K^y 



d^z 



dT,zz 



(82) 



d{qlql) 



247r2, 



dqt 



These expressions reduce to those found for the rigid the- 
ory when V = 0. 

The calculation of one-loop diagrammatic corrections 
to Tyz and F^j^ is similarly straightforward. Ey^ is given 
by the diagram in Fig. ^(a). This amplitude is propor- 
tional to V since it is formed by contracting dyUy{dxUz)^ 
with dzUz{dxUz)^- ^yy is given by the diagram in 
Fig. ^(a) ; it must be proportional to tJ^ since it is formed 
by contracting dyUy{dxUz)^ with itself. The one-loop cor- 
rections to Tyz and Tyy are given below to lowest order 
in e: 



diqyQz) 



gv 



1 



(83) 



yy 



dql 



=0 



We then use the conditions imposed on the vertex func- 
tions in (HI) and (^l|) and the one-loop diagrammatic 
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corrections in ( |82D and (|8^) to find the renormalization 
constants (the Z's) in terms of g and v. We find that the 
relations for Zy, and Zg are unchanged to zeroth or- 
der in V. Z and Zt^ also have terms that are independent 
of V as shown below to lowest order in e: 



Z 

Z., 



1 



1 



9 



127r2 



(84) 



The variation of g and v with the length scale /i is ob- 
tained by enforcing that both bare coupling constants 
do not depend on /z, i.e., we set ^dw/dfi = ^dv/dfi = 0. 
These two requirements determine the recursion relations 
for g and v; we find that dg/dl is unchanged to lowest 
order in v and 



gv 

167r2 



(85) 



The zeroth order solution for g was found previously in 
(|6^); we plug this solution into (^5|) and find 



v{l) = 



[l + 5oV(67r2)] 



3/8 ' 



(86) 



where vq = Byy^/BvB^ and go = JB^/K^y/Kl 



C. Renormalized Elastic Constants 

We found in the previous two sections that the renor- 
malized elastic constants are obtained by solving the 
Callan-Symanzik equation for the renormalized vertex 
function. We find the CS equations for using the fol- 
lowing scaling equations which relate the bare and renor- 
malized vertex functions: 

r,,(q,u;,i;) =Z-i/32^rL(q',g,tJ,^) (87a) 
r,,(q, «;,!;) = Z'^ Z'^ Z^'^Tly{c{ , g,v, ^i) (87b) 

VyM.w.v) = z-i/2z^r;,(q',g,ij,^). (87c) 

Eq. ( ^7a| ) yields a CS equation identical to (^) to lowest 
order in v, and thus the renormalized elastic constants 
i3^(q), K^^{(\}, and K^y{q) are identical to those ob- 
tained in ( |67| ) using the Uy = theory. The fact that 
the elastic constants are identical to zeroth order in v is 
a consequence of the fact that the nonlinear term pro- 
portional to TJ does not introduce any harmonic terms 
that were not already present in the theory without Uy 
fluctuations. We also find that the coefficient of T'[^y{c^) 
is unity to lowest order tJ, and hence the vertex function 
Vyy does not rescale. As a result, — By {I — 0) plus 
higher order terms in v. 

We do, however, find a n ontri vial renormalization of 
S^^. The scaling relation in (87c) leads to a CS equation 
for with a similar form to the one found in (|6^) . We 
find 



via) _ vigl f _d_ 

2 3 \'^'dq. 



= 



to zeroth order in v, where 

d(lnZ) 



dg 



127r2 



(89) 



and 7/ and rjy were defined previously. The solution to 
( ^ ) can be transcribed from ( |6l| ) and is displayed below: 



rL(q,5,i'(ff),/^) = exp 



rjydl' 



(7j^,exp 



3./0 



dl' 



T]dl' 



(90) 



qz,9,^j.o 



Since 77, rjy, and r] scale as l/l, the integrals in the ar- 
guments of the exponentials scale logarithmically with 
I. Thus, the exponentials yield power-laws in g, and we 
find, for example. 



exp 



dl' 



9{l) 



1 5/8 



90 



(91) 



The renormalized vertex function in ( |90| ) obeys a scaling 
form analogous to the one obeyed by the renormalized 
sliding columnar vertex function in (p4). We find 



rL(q,5,Ai) = b ^ [g/go 



5/8 



(92) 



'^Iz (bqx,bqy [g/go] ^h^q^ [s/so]^^*^ , 5, A^o' 



where the 6 ^ prefactor is present because y scales as 
h and z scales as 6^. We then choose h = fx^^ = [q^ + 
qlql+w~^q^]~^^'^ = [/i(q)]~^ to match the conventions of 
the previous section, substitute ( |6^ ) for g/go, and return 
to variables with dimension. The renormalized vertex 
function becomes 



2By^ 


\ , -90 1 















-3/4 



qyQz, (93) 



where /I and /i(q) were defined previously. The renor- 
malized elastic constant i?*'^(q) is the coefficient of qyq^ 
in the above expression. Therefore, we find that both 
B^ and B^^ scale to zero logarithmically with q at long 
wavelengths defined by h{q) ^ A^^"^ expl—Gn^ / go]. 



V. CONCLUSION 

We have calculated the Grinstcin-Pelcovits renormal- 
ization of the elastic constants for the sliding columnar 
phase. We first used a simplified model of the sliding 
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columnar phase in which the DNA columns were pre- 
vented from fluctuating perpendicular to the lipid layers. 
We found that the elastic constants scaled as powers of 
ln[l/g] at long wavelengths. In particular, we found that 
the compression modulus B scales to zero and the rota- 
tion and bending moduli Ky and K scale to infinity as 
q tends to zero. We then added in perpendicular fluctu- 
ations of the columns perturbatively and found that the 
above results were unchanged to lowest order in the cou- 
pling between strains parallel and perpendicular to the 
lipid layers. We employed dimensional rcgularization in 
our RG analysis of the sliding columnar phase to ensure 
rotational invariance. RG schemes that break rotational 
invariance, such as the momentum-shell technique, did 
not yield correct results. 
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(a) 




(b) 



u 




u 



FIG. 5. Schematic representation of the two relevant non- 
linear terms in both the 3D smectic and sliding columnar 
elasticity theories. The perpendicular derivatives (±) corre- 
spond to the 3D smectic theory and the x derivatives to the 
sliding columnar theory. The term (c^x,^,!!)'* is pictured in (a) 
and the term {dzu)(dA_,xuf is pictured in (b). The symbols 
_L, X, and z represent _L, x, and z derivatives of the u field. 
The diagram with four u fields in (a) does not contribute to 
the renormalization to one-loop order; only contractions of 
(b) with itself contribute. 



APPENDIX A: EVALUATION OF THE 3D 
SMECTIC ONE-LOOP DIAGRAMS 

Our t ask in this Appendix is to calculate S(q) deflned 
in Sec. [I C as the one-loop diagrammatic corrections 
to r(q), the vertex function for the 3D smectic. These 
corrections arise from the nonlinear terms in the Hamilto- 
nian in (]Io|). The two nonlinear terms are dzU^CV /2 
and (Vxm)*/8 (shown schematically in Fig. |[), and only 
contractions of the former contribute to the renormaliza- 
tion to one-loop order. The three possible contractions 
are shown in Fig. |^. The diagrammatic corrections S(q) 
can be expressed as 

I](q) = ni(q)g2 + Halq)?! = Ei(q) + S2(q). (Al) 



Note that we have separated the and 
of S(q) so that to lowest order in q 



dependence 



and 



'dqi 



d^2 
dqi 



(A2) 



The contributions oidY^-z/dql to dY./dq'^ and of rfSi/dq^ 
to dS/dg^ at the special point q^ ~ and q^ = are 
higher order in e than the contributions in ( [A^ ) and ([A^). 
We begin by calculating Si(q). 




(C) 



-L, jc 




-L, X 



-L , X 



FIG. 6. The three one-loop diagrams that contribute to 
the renormalization of the 3D smectic and sliding columnar 
elastic constants. These diagrams are formed by contracting 
dzu{d±^xu)^ with itself. The diagram in (a) contributes terms 
proportional to since a factor of is on each external leg. 
The diagrams in (b) and (c) contribute terms proportional to 
(A3) <1± in the 3D smectic theory and terms proportional to q^Qy 
and in the sliding columnar theory since these diagrams 
have q^ or q^ on the external legs. 



1. Calculation of Si(q) 

The diagram in Fig. ^(a) alone contributes to Si(q) 
since it is the only one with on the external legs. To 
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evaluate the integrals in the perturbation theory, we use 
dimensional regularization, i.e we take d = 3 — e, set 
the cutoff to infinity, and look for the 1/e terms. Si(q) 
is obtained by calculating the contribution fi'om the 
following integral: 



Si(q) 



9£ 
2 



(p-'k r 



-oo (27r) 



3-e 



fci,fci,G(k + q)G(-k) 
where i, j = x,y and 
G(q) 



{q± + kj_)i{q± + k±)j X 

(A4) 



1 



— 1 4 ' 



(A5) 



The coefficient of the term in (A4) is ni(q). We can 



then approximate Si(q) by writing Si(q) = q^Ili{qj_ = 
0, qz) plus higher order terms in tha t va nish when we 
apply the boundary condition in Eq. (23a). We obtain 
^i{Qz) by s etting q± ^ in the integral on the right hand 
side of ( [A4| ). 

To evaluate the integral, we first combine the denom- 
inators of G(k + q) and G(— k) employing the following 
identity: 

1 1 



dx- 



kl 
1 



{kz + xqzY x{l — x)q1 + w ^fc^ 



(A6) 



We then change variables to k'^ = k^ + xqz and perform 
the integration over fc^. We find that Si(q) can be writ- 
ten in terms of the integral J(4, 3, x, qz) with J(s, v, x, qz) 
defined by 



J{s,v,x,qz) = 



dk^k]^"' 



, r -(2v 
AT{v 2) V4 



x{l — x)q'^_ + w ifc^]"^^ 
-2 + e) 



xr( -Gs + 2-e) 



X [x{l ~ x)wq 



2-1 (s-2ti+2-e)/4 



(A7) 



where T{x) is the gamma function evaluated at x. The 
expression for Si(q) is simple when expressed in terms 
of the integral J(4, 3, x, qz); we find 



Si(q) = / dx J{4,3,x,q, 

167r 7o 



(A8) 



From (A7) we know that the most dominant term in 
J{4:,3,x,qz) scales as 1/e and thus 



Si(q) = - 



{ml) 



-e/4 



plus higher order terms in e. We can also write 5]i(q) as 



dSi(q) 



dql 



(AlO) 



when we replace w by (g/i*^) 



2/3 



2. Calculation of S2(q) 

E2(q) is determined by calculating the contribu- 
tions from the diagrams in Figs. ^ (b) and (c). S2(q) is 
the g^ part of the the following integral: 

/d^-'^k r 
{kz + qz){k± + q±)jkj_,kz 



G(k + q)G(-k). 



(All) 



The contributions come from expanding G(k -I- q) to 
second order in q±; we see from (All) that we need 
both the first and second order terms in the expansion. 
The coefficient of the term in the above expansion is 
n2(g± = 0,qz), and thus S2(q) = q\ll2{qz) plus higher 
order terms in q± that vanish when we apply the bound- 
ary condition in ( ^3b| ). 



The first and second terms in the integrand of (All) 
correspond to the diagrams in Figs. ^ (b) and (c), re- 
spectively. We break up the integral so that S2(q) = 
^2(1) + ^2(1)1 ^■'^d we first calculate S2(q). 



S^(q) = - 



1 

dkz f di} 
2^ 



(A12) 



(2^) 



2-£ 



dk I k 



{kz + qz)^k^,k^jG{-k) 



d^G{k + q) 



dq^idqi 



where VL is the solid angle in 2 — e dimensions and the sec- 
ond derivative of G gives the coefficient of the quadratic 
term in the expansion of G(k + q). We then remove the 
angular dependence by integrating over f2 and using the 
following two identities: 



fc I ,-fc I = — ^— ^fc^ 5. 



and 



(27r)2-^ ^'^^ 2-. 



I. I. i. I. 



(A13) 



(27r 



(A14) 



5*2 

rk^{dijSi 



(2-e)3 

where Sij is the Kronecker delta and Sd — il/(27r) 



(A9) 27r'*/V((27r)''r(rf/2)) with d = 2 - e. We are interested 
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in the lowest order terms in e and hence will use 5*2-6 ~ 
(27r)~^ below. We then change variables to k'^ — + Qz 
and combine the denominators of G(— k) and G(k + q) 
using an identity similar to (A6). 



1 



1 



{kz-qzy + w-^ki 

»l 



fn{x) 



(A15) 



r{n+l) / dx 

■^0 \{kz-xq,)^ +x{l-x)ql+w-'^kl 

where n = 2, 3 and 



n+l ' 



fn{x) = 



1 



(l-a;)V2, n = 3. 



(A16) 



We change variables again to fc" = A;^ + x(7^ and integrate 
over fc"; we find that S^Iq) can be written in term s of 
the integrals J{s,v,x,qz) defined previously in (A7): 



327r 



dx 



-5(l-x)J(4,3,a;,(7,) 

-l?>x^{l - x)qlJ{4:,b,x,qz) 
+9w;-\l - 2:)V(8, 5, 9^) 

+45^-^x2(1 - xfqlJ{S, 7, X, . (A17) 



J(4, 3, X, qz) and J(8, 5, qz) have terms proportional to 
1/e but J(4, 5,x, g^) and J{8,7,x,qz) do not. We keep 
the terms that are proportional to 1/e and drop the oth- 
ers. In the last step we perform the x integration and 
find 



,,,1/2 



2N-e/4 



(A18) 



plus higher order terms in e. 

We next obtain '^2(1) by calculating the contribu- 
tions from the diagram in Fig. ^(c) . S2 (q) can be written 
in terms of the following integral: 



dG(k + q) 



(2^) 

kz{kz + qz)G{-'k) kj_,q±jq^^ 



dqi 



9_L=0 



q±,q±^ d^Gik + q) 



dqj_^dqA 



q±=o 



(A19) 



The first and second derivatives of G give the coefficients 
of the linear and quadratic terms in q± in the expansion 
of G(k -I- q). We then follow a procedure similar to the 
one employed to find Yl'^i'i)! '^6 change variables to 
k'^ = kz + qz, combine the denominators of G(k + q) and 
G(~k ), and integrate over O. The remaining integrals in 
( A19| ) are over k± and x. We then integrate over k± and 



write S2(q) in terms of J(s, v, x, qz); we find 



—qi / dx -9(l-a;)J(4,3,a;,g.) 
327r Jq L 



+27x{l-xfqlj{4,5,x,qz) 
+9w-^{l~ xfj{8,5,x,qz) 

-45w"^a;(l - x)^qlJ{8, 7, x, qz) 



(A20) 



Only J{A,3,x,qz) and J(8,5,a;, g^) have terms propor- 
tional to 1/e. We keep these terms and perform the in- 
tegration over X to find 



(A21) 



We obtain E2(q) by adding E^(q) and S^(q) in ( |aT^ ) 
and ([A21D to yield 



dS2(q) 



dq 



-2/3 g 

327re' 



(A22) 



once we set w = {g^'^)'^^^ and ignore higher order terms 
in e. 



APPENDIX B: EVALUATION OF THE SLIDING 
COLUMNAR LOOP DIAGRAMS 

The aim of this Appendix is to calculate S(q), the 
one-loop diagrammatic corrections to the vertex function 
for the sliding columnar phase. The rotationally invari- 
ant theory given in ( p5| ) contains two relevant nonlinear 
terms, dzUzidxUz)"^ and {dxUz)'^- These terms are pic- 
tured schematically in Fig. |5[ From this figure we see 
that only contractions of dzUz{dxUz)^ renormalize the 
elastic constants to one-loop order. The three possible 
contractions are shown in Fig. ^. S(q) has ql, q^qy, and 
q^ contributions, and we will calculate each separately 
below. To do this, we express S(q) as 



E(q) - ni(q)g2 + n2(q)g2g2 + n3(q)g^ 
= S]i(q) + S]2(q)-f I]3(q). 



(Bl) 



We have separated the q^, q^Qy, and q^ dependences so 
that, for instance. 



d4 



(B2) 



As in Appendix we use dimensional regularization to 
calculate the integrals. 



1. Calculation of Si(q) 

The contribution to S(q) results from squaring the 
diagram pictured in Fig. ^b) and contracting both pairs 
of X derivatives. This leaves qz on each external leg as 
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shown in Fig. ^(a). Si(q) is the q1 part of the following 
integral: 



Sl(q) 



where 



9£ 
2 



(27r) 



3-e 



(g. + fcx)'fc^G(q + k)G(-k) 



G(q) 



1 



(B3) 



(B4) 



The coefRcient of the in the above integral is Hi (q) and 
thus Ei(q) = q^\ii{qx,y = 0, plus higher order terms 
in qx and qy that vanish when we apply the boundary 
condition in (|55|). Thus, Si(q) is obtained by setting 
Ix = % = ^ in (|B3|). We find 



Sl(q) 



w 



-1/2 



9^ 



2 (27r)3-^ 

1 



(B5) 



-11,2 1,2 

ri/,„ ri I 



w 



-11,2 1,2 



w ^^^fcy and 



where we have changed variables to ky 
dropped the prime. The first step in evaluati ng t his in 
tegral is to combine the two denominators in (|Bq) using 



the identity in (A6) with fc^ replaced by k'^k\. Wc then 



perform the integration over kz and find that Si(q) can 
be written in terms of the integral /(4, 0, 3, x, q^), where 



I{s,t,v,x,qz) 



dkx dky 



l,s Lt — e 
l^x'^y 



x{l ~ x)q1 + w ^k'^k\ 



v/2 



8r v/2 yr 



T{-{s-t + e) 



x{l — x)wq^ 



t - s - 2 + e) 

(s+t-2u+2-e)/4 



2. Calculation of E2(q) 

Both the q^qy and gj^ contributions to S(q) come from 
the diagrams with x derivatives on the external legs. The 
two contributing diagrams are shown in Figs. |^ (b) and 
(c) . Their sum is given by 



Sum = 



(BIO) 



(f-^k 



(fc, + qzYki + {qz + kz){qx + kx)h 



(27r)3 
xG(k + q)G(-k) 



We find the q^qy terms by expanding G(k + q) to second 
order in qy. We see that only the quadratic term in the 
expansion contributes. Higher order terms will vanish 
when we apply the second boundary condition in js^). 
We then follow a procedure similar to the one employed 
to calculate the contribution to the 3D smectic vertex 
function in Appendix Wc find that S2(q) can be 
written in terms of the integrals I{s,t,v,x,qz) as shown 
below: 



S2(q) 



Stt 



1/2 

2 2 

2 ^x^y 



dx 



-2(l-a;)/(4,0,3,a;,g,) + 



6w-\l-xfl{6,2,5,x,qz)~ 
ixql{2x - 1)(1 - a;)/(4, 0, 5, x, q^) + 

l^w-^xql{2x - - xfl{Q,2J,x,qz)\. (Bll) 



We look for the leading order terms in e in (Bll); 
/(4, 0, 3, a;, gz) and /(6, 2, 5, (?z) have leading order 
terms proportional to 1/e while /(4, 0, 5, x, gz) and 
7(6, 2, 7, X , qz) do not and are dropped. After integrat- 
ing (Bll) over x we obtain 



and 



(B6) 



S2(q) = 



247r26 



<Ix%{mz 



2N,-e/4 



d{qlql) 
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247r2e ' 



(B12) 



(B13) 



We give the most general form for the integrals over kx 
and ky since we will need these integrals later when we 
calculate S2(q) and S3(q). We find 



Si(q) - 



-7/,-l/2 



87r2 



-qz I dx 1(4,0,3, x,qz). 
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and 



Sl(q) 



8^2e 



qliml)-'" 



(B7) 



(B8) 



since 7(4, 0, 3, gz) cx 1/e. We then set w — gji'^ to find 
I]i(q) as a function of g, 

dSi(q) 



3. Calculation of S3(q) 

E3(q) is obt ained by calculating the terms propor- 
tional to q^ in ( p31Cl| ). We obtain these terms by expand- 
ing G(k -I- q) to second order in qx and noting that both 
first and second order terms in the expansion contribute. 
Note that higher order terms in the expansion will van- 
ish once we apply the third boundary condition in ([s^). 
We calculate the g^ contributions from Figs. || (b) and 
(c) separately and define S3(q) = £3(0) -I- S3(q). We 
(B9) first calculate the contribution from Fig. ^(b). Using the 
same procedure as the one employed to calculate the g^g^ 
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contribution to E(q), we find that S|(q) can be written 
in terms of the integral I{s,t,v, x, Qz). 



,,,-3/2 



S§(q) 



, dx - (l-x) (6/(4, 0,3, a:, 9,) + 

/(2, 2, 3, X, q,) + l^x\ll{^, 0, 5, X, q,) + 
Sx\ll{2,2,5,x,qz)) + 

3w-^{l - xf (4/(8, 0, 5, x, qz) + 

20a;2g2/(8,0,7,2;,g,) + 

4/(6, 2, 5, x, q,) + 20x2q2/(6, 2, 7, x, q,) + 

/(4, 4, 5, x,qz) + 5x\ll{4, 4, 7, x, g,)) ] . (B14) 



We note that three of the integrals in (B14), 



/(4,0, 5,x, g^), /(8,0, 7,a:, g^), and /(6, 2, 7, x, g^), have 
leading order terms that scale as and are dropped. 
Two integrals, /(2, 2, 3, x, g^) and /(4, 4, 5, x, gz), have 
as well as 1/e terms, while the remaining five 
integrals /(4, 0, 3, a;, g^), /(2, 2, 5, x, g^), /(8, 0, 5, x, g^), 
/(6, 2, 5, X, qz), and /(4, 4, 7, x, g^) have leading order con- 
tributions that scale as 1/e. We collect terms and per- 
form the X integration to find 



1 



^'M = --^qt{wql)-^^^ 



1 + ln[2l - — 
e ^ ^ 12 



(B15) 



Note that the dominant contribution to Eg is of order 
e"^ rather than e~^. The undesirable term and 
the ln[2]/e term will be cancelled by terms in E§. The 
term proportional to ln[2]/e originates from the integrals 
/(2, 2, 3, X, g^) and /(4, 4, 5, x, g^). This can be seen by 
expanding /(4, 4, 5, x, g^) in powers of e; we find 



/(4,4,5,x,g,) 



2u;-^/2 



[x(l — x)wq 



2 r(5/2) 

21-^/4 



2 r(i) 



to order 0(l/e), where r'(x) is the derivative of the 
gamma function evaluated at x. The logarithm arises 
from evaluating the derivative of the gamma function at 
a half integer. For example, r'(5/2)/r(5/2) = -7+8/3- 
21n[2] where 7 is the Euler-Mascheroni constant. 

We can also write E3(q) in terms of the integrals 
I{s,t,w,x,qz). We obtain 



-3/2 



2 



8^^ .0 

,2 i 



dx 



(1 - x)( -10/(4,0, 3, x,g,)-h 

30x(l - x)g^/(4,0,5,x,g2) - 
3/(2, 2, 3, X, qz) + 9x(l - x)gf /(2, 2, 5, x, g,)) + 

3u;"^(l - x)2 (^4/(8, 0, 5, x, qz) - 

20x(l - x)g2/(8, 0, 7, x, g,) + 4/(6, 2, 5, x, g,) - 
20x(l - x)g2/(6, 2, 7, x, g,) -t- /(4, 4, 5, x, g,) - 

5x(l-x)g2/(4,4,7,x,g,))], (B17) 



which becomes 
^3(^) = -8;J^ 



-QxiwQz 



2N-6/4 



1 7 

ln[2l 

e ^ ^ 12 



(B18) 



when only terms proportional to 1/efand 1 /e a re re- 
tained. We see that when we add ( |B15| ) to (|B18| ), the 
terms proportional to 1/e^ and ln[2]/e cancel and we are 
left with 



S3(q) 



127r2, 



-.lii^q^z 



-£/4 



and 



dS3(q) 



dqt 



9 



(B19) 



(B20) 



APPENDIX C: FINITE WAVENUMBER CUTOFF 

In this Appendix we show that employing a finite cutoff 
leads to ambiguities when we evaluate the sliding colum- 
nar one-loop diagrams. These diagrams are shown in 
Fig. H; (a) contributes to Si(q) and both (b) and (c) 
contribute to T,2{q) and S3(q). The ambiguous result is 
that we obtain different answers for S(q) depending on 
whether external momentum q is sent through the top or 
bottom part of the internal loop. The ambiguity devel- 
ops when momentum q^ appears in the internal loop and 
the top and bottom paths through the internal loop are 
different. The diagram that causes this ambiguity is the 
q^ part of Fig. ||(b). We can see this by calculating the q* 
corrections to the vertex function, E|(top) and E§(bot), 
which result from sending k + q through the top(bottom) 
sections of the internal loop. 



(B16) E^(top) = 



d^k 
(2^ 



[fc'(fc. + g.)'G(-k)G(k + q) 



(CI) 



and 



S^(bot) 



d^k 

[kl{kz+qz)'G{-k)G{k + <i)], 



(C2) 

where A is a finite wavenumber cutoff and G{q) was de- 
fined previously in (|B4[). With A ^ 00, 



(C3) 



E^(top) ^ S§(bot). 



If we employ dimensional regularization instead and send 
A ^ CX3, these top and bottom amplitudes are identical. 
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APPENDIX D: DERIVATION OF THE 
NONLINEAR STRAINS IN THE PRESENCE OF 
FLEXIBLE MEMBRANES 

In this appendix, we derive expressions for the non- 
linear strains Uyy{'x.) and Uz^(x) introduced in (70a) 
and ( 70b ) for the case of flexible membranes. A com- 
plete description of lamellar DNA-lipid complexes re- 
quires separate coordinates for each membrane and each 
DNA molecule. Displacements of membranes and DNA 
molecules parallel to the membrane normals (along the 
y-direction when the membranes are flat) are locked to- 
gether. We can, therefore, model the complexes as a 
stack of membranes each with a one-dimensional mass- 
density wave representing the DNA lattice just above 
it. We employ mixed Lagrangian-Eulerian variables in 
which the coordinate y = na specifying the layer or mem- 
brane number is a Lagrangian variable and the coordi- 
nates (x, z) = r are Eulerian variables specifying posi- 
tions in a flxed projection plane. The positions of mass 
points on membrane n are then given by 

R„(r) — XX + zz + [na + Uy{na, r)] y. (Dl) 

The density in membrane n can be expanded as Pn(r) = 
Pn + V'n(r) + V'n(r), whcrc is a constant, '0n(r) — 
|V'„|e*'^"('"), and 



0„(r) ^ko[z- u^{na,r)] 



(D2) 



with ko = 2TT/d. 

To construct the strain variable Uyy{x.) with x = (y, r), 
we introduce the distance /„(r, r') between points r on 
membrane n and r' on membrane n + 1 via 



/2(r,r') = |R„+i(r')-R„(r)r 



The shortest distance between a point r on membrane n 
and any point on membrane n -|- 1 is then 



Z2(r)=min/2(r,r'). 



(D4) 



where = • is the metric tensor of the nth surface 
and where we used y ■ = d^Uy. We then minimize 
Z^(r,r') over 5r^^ and obtain 



a{l + dyUy) g^^di,Uy 



(D8) 



and 



^2 , (r) = (1 + Q^UyY (1 - g^'^d^uyd^uy) . (D9) 



Finally, using g^" — (5^1/) ^ where 



we obtain 



Uyy{x) — 



(l+dyUy)^ 

1 + {Vuyy 

1 



(DIG) 



(Dll) 



~ 9yUy ~ - [(d^.Uy)'^ + {d.Uyf - (dyUy)'^] , 

with V = {dx, 0, dz). It is straightforward to verify that 
Uyy{x.) — for a uniform rotation of the entire system. 
For example, a rotation of the system about the z axis by 
9 produces strains dyUy — 1/ cos9 — 1 and d^Uy = tan 6' 
which to vanish. 

The strain Uzz (x) can also be defined in a rotationally 
invariant way via 



.(x) = ^ [fco'-5'"'5^'/'(x)9.0(x) 



(D12) 



where 0(x) 

order in d^Uz and d^^Uy, the nonlinear strain Uzz is 



y/a{^) is defined in (D2). To quadratic 



(D3) Uzz{y^) w dzUz 



{dxUzf + {dzUzf - {dzUyf , (D13) 



where ^^(x) = u^(y,r). 



The strain variable Uyy is defined as 



^yy 



(D5) 



This quantity is by construction invariant with respect to 
global rotations of the entire system. To evaluate Uyy (x) , 
we expand R„+i(r') — R„(r) to lowest order in 5v ~ r' ~v 
and a: 

R„+i(r') = R„(r) + a (1 + dyUy{^)) y + Sr^e^, (D6) 

where fi — x,z, = 9pR„(x) is a covariant tangent- 
plane vector of the nth surface, and Uy{x.) — Uy{na,r). 
Then 

ll{r, r') = (1 -f dyUyf 2a (1 + dyUy) 5r''d^Uy 

+ g^.Jr^^5r\ (D7) 
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